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Abstract 



A consequence of non-Gaussian perturbations on the Sachs- Wolfe effect is studied. For a 
particular power spectrum, predicted Sachs- Wolfe effects are calculated for two cases: Gaus- 
sian (random phase) configuration, and a specific kind of non-Gaussian configuration. We 
obtain a result that the Sachs- Wolfe effect for the latter case is smaller when each temper- 
ature fluctuation is properly normalized with respect to the corresponding mass fluctuation 
tjf(R). The physical explanation and the generality of the result are discussed. 



I. INTRODUCTION 

It is a challenging problem for cosmologists to explain the large scale structures. It is gen- 
erally believed that some microphysical processes produced small amplitude perturbations in 
the early universe, and that the present large scale structures grew by gravitational instabil- 
ity of these perturbations. Perturbations are called Gaussian (G) perturbations if the phases 
of each fourier mode are random, and they are called non-Gaussian (NG) perturbations if 
otherwise. 

G initial fluctuations are well motivated in the context of inflationary cosmology. During 
the inflationary phase, quantum fluctuations of a scalar field, which drives the inflation, 
produces G fluctuations [1-4]. There have been many, both analytical and numerical, studies 
of evolution of G fluctuations [5-7]. An evolutionary scenario requires further assumptions 
about the type of matter in the universe, (e.g., Cold Dard Matter (CDM) or Hot Dark 
Matter (HDM)), and the initial power spectrum of the fluctuations p(k) = \5k\ 2 - In some 
cases even fl = pj ' p critical and A, the cosmological constant, are left as free parameters. 

Statistical methods are typically used to make comparisons between predictions of a 
model and the observations. Here a further assumption about the relationship between the 
distribution of matter and the distribution of observed galaxies needs to be made [8-10], i.e., 
biased galaxy formation. 

For a model to be branded a success it must fulfill at least two requirements. First, it 
must be capable of reproducing the observed large scale structures. Second, its predicted 
Cosmic Microwave Background Radiation (CMBR) must agree with the observations [11-14]. 
There are several mechanisms contributing to the CMBR temperature fluctuations, but for 
fluctuations on large angular scales {6 ^> 1°) only two are important. The first is the doppler 
effect coming from our pecular velocity with respect to the CMBR rest frame. The second 
is the Sachs- Wolfe effect, which is caused by fluctuations in gravitational potentials on the 
last scattering surface [15]. With the recent results from the Cosmic Background Explorer 
(COBE) [14], some G models are on the process of being ruled out. For example, HDM, 



Q — 1, with the Harrison-Zel'dovich power spectrum predicts too large =jj- in direct conflict 
with COBE [16,17]. 

An alternative to Gaussian models are whole class of non-Gaussian (NG) models. Exam- 
ples of NG models are Cosmic Strings, Monopoles, and Textures, which are the remnants of 
phase transitions at early universe [18-23] . Production of some NG fluctuations can even be 
accomodated within the context of inflationary scenario [24-26] . There are also mechanisms 
for generating NG perturbations without resorting to 'exotic' physics of the early universe, 
for example the explosion scenario of Ostriker and Cowie [27], and Ikeuchi [28]. 

Detailed N body studies for evolution of NG perturbations have been done [29-33] , and 
their corresponding Sachs- Wolfe effects were calculated [34,35]. In this paper, we would like 
to do a comparative study of Gaussian vs. non-Gaussian fluctuations in the context of Sachs- 
Wolfe effect. The question we would like to address is whether the predicted Sachs- Wolfe 
effect can be reduced by describing the primordial density perturbations by a particular NG 
field configuration instead of a G field configuration. That is, for a fixed power spectrum 
p(k) = \Sk\ 2 , we would like to consider the consequences of two types of field configurations, 
^0*0 = J2 \5k\e t ^e lk ' x . The first is the G case where all the phases e 1 ^ are random. The 
second is a NG case where the phases are correlated in some way. For each case we will 
calculate the predicted ^f- due to the Sachs- Wolfe effect and compare the magnitudes after 
normalizing each ^f- by the corresponding mass fluctuation ^ (R) at some scale R. 

Even for a fixed power spectrum it is obvious that the set of all possible NG field con- 
figurations is an infinite dimensional parameter space. There have been some analytical 
work on the parametrization of NG perturbations ( Moscardini et al. [36] , Weinberg et al. 
[37]). Because their parametrization is in the 1 point probability function space and not 
in fourier or configuration space their method is not directly applicable here. The reason 
is that every possible field configuration described by a different set of phases e 1 ^ gives 
a probability function P[S], but each probability function does not give a unique field con- 
figuration. Therefore, in this paper, we will first study the consequences of a very special 



NG field configuration, i.e., e 1 ^ are all +1 (or all -1). Clearly the resulting configuration 
— (x) is a spherical overdensity (spherical void). The difference between this NG and the G 
configuration is most accute at the center of this object, where all the amplitudes of different 
wavelength fluctuations add linearly in the NG case. 

In this paper, for this NG field configuration, we have done a model calculation for the 
power spectrum p(k) ~ 5{k — k Q ). As will be shown in detail, we arrive at the result that the 
predicted Sachs- Wolfe effect for this NG case, when properly normalized by mass fluctuation, 
is smaller than the corresponding prediction for the G case. 

The layout of the paper is as follows. In section 2, some notations and basic formuli are 
quoted from literature. In section 3, and ^jf(R) are calculated for the power spectrum 
p{k) ~ 5{k — k ) as G fluctuations. In sections 4 and 5, for the same power spectrum, we 
calculate ^ and 5 -§(R) again, respectively, but now as NG fluctuations. In section 6, the 
calculated temperature fluctuations for the two cases are normalized with respect to the 
corresponding mass fluctuation, and they are compared. In section 7, physical explanation 
of the result is given and some comments and generalization are given. 

In the model calculations the scale factor of the universe and the density perturbations 
were assumed to evolve as in flat CDM universe. 

II. PRELIMINARIES 

By definition, a spatial perturbation is just a fourier sum of fluctuations on all wave- 
lengths. 

5 JL{g) = ]T |%|e^ fc V fo (2.1) 
P k 

A power spectrum is defined as 

~~ * 2 ( ~^ J ^\ 2 "-i 

p{k) = \5j:\ — ► I — J \5%\ d k continuum limit. (2.2) 

L is an arbitrary length which is large compare to any length scale in the problem. If the 
fluctuations were produced by some microphysical processes in the early universe, then, a 



priori, there should be no preferred direction. Therefore, rotational invariance dictates that 

— * 

p(k) =p(\k\). 

A model power spectrum we will use is 

P (k) = \ 6 ° ^ 2 I - fc - + 2 I (2.3) 

y otherwise. 

We will not use the Dirac delta function, i.e., p(k) ~ 8(k — k Q ), because for a NG case we 
will need to evaluate integrals of the form / \5 k \d 3 k{...}, and square root of a delta function 
is tricky to handle. We will assume that AA; is thin enough so that / f(k)dk ~ f(k )Ak for 
any reasonable function f(k). 

It turns out that this is the simplest power spectrum with which we can demonstrate the 
essential differences between the G and a particular NG configuration, which we will next 
discuss. 

With this power spectrum, we will consider two different field configurations, i.e., choice 
of phases. 

random, Gaussian case (G) , , 

all + 1 (or all — 1), a non — Gaussian case (NG) 

The physics behind generating a CMBR fluctuation from a density fluctuation is easy to 
understand. Simply stated, a matter perturbation distorts space, and therefore it 'distorts' 
geodesies of photons which were last scattered at recombination. Because a Newtonian 
potential is time independent in a linear regime only the perturbation in the Newtonian 
potential on the last scattering surface is important in the linear Sachs- Wolfe calculation. 
Contributions from any intermediate distortion of a geodesic is unimportant by the blue and 
redshift cancellation effect. If a structure becomes nonlinear, then the gravitational potential 
is no longer constant, and the cancellation is not complete. This higher order effect has been 
studied by Rees and Sciama [38]. 

The relationship between a matter perturbation and ^ can be described most succintly 
by the expression [15, 39] 




= m | T ° = I mx+r q) - sm] , (2.5) 

where S(f>(x) is the Newtonian potential at point x, r Q = 2if~ 1 , q is a direction in sky, and 
T Q = 2.7K is the mean CMBR temperature. From V 2 <50 = 4nG5p, we get 

50(f) = -\Hl j: 8 k {t )k-*<&*. (2.6) 

Rigorously, these expressions are valid in time-orthogonal coordinate and for a flat universe. 
The generalization to an open or closed universe is discussed by Anila and Motta [40]. 

III. GAUSSIAN FLUCTUATIONS 

For G fluctuations, because of the random phase of each 5^ (2.1), only statistical infor- 
mation can be studied, e.g., the angular correlation function of Using (2.5) and (2.6) 
we get 



,AT AT 
— (q u x) — 

L\ 3 rdk 



C(qi-q 2 ) = (—(qi,x) — (q 2 ,x)) space 



(Li \ ( dk 
2tt ) J ~¥ V ^ tioikTolq! - q 2 \) - monopole - dipole) . (3.1) 

The monopole (i.e., spherically symmetric) term is subtracted out because it is unobserv- 
able for a fixed observer. The dipole term is also subtracted out because in principal it is 
impossible to distinguish an intrinsic dipole anisotropy from the observer's pecular velocity. 
The monopole and the dipole terms can be identified by using the following identity, which 
is proven in the appendix. 

U\AA - BB\) = £(2Z + l) 3l (A)UB) Pl (A ■ B) (3.2) 
i 

j and pi are Spherical Bessel and Legendre functions, respectively. Setting A = kr Q , A — q ± , 
B = kr , and B = q 2 we identify the monopole and the dipole terms to be jl{kT ) and 
3ji(fcT ) 2 pi(<?i • q 2 ), respectively. 



Using a thin Ak approximation to evaluate the integral for the power spectrum (2.3), we 
finally get [15,41] 



C(& • ft) = nHt(£f^ (^-fo(koT ) - mkor fpMi ■ &)) , (3-3) 

with c = r |ft — ft|. The overall normalization <5 D must be fixed. This is usually done 
by relating it to either ^jf(R) or mass autocorrelation function £(r) = (x + r)^{x)) space . 
Normalizing <5 G by £(f) would be subjecting the model to more stringent requirements, which 
would be appropriate for a more realistic power spectrum like p(k) ~ k n . Therefore, for our 
simplistic models, we choose to normalize it by ?jf(R). 

Again for the G case, we will just quote the well known result relating the power spectrum 
and *£{R) [39]. 

(\ 5 -^(R)\ 2 ) s P ace = (^)7 d*kp(k)W(kR) (3.4) 
where W{kR) is a window function. We will use the top down window function 

W(y) = 9y~ 6 (siny — y cosy) 2 . (3.5) 
Evaluating the expression for the power spectrum (2.3) we get 

(\ 6 -^(R)\ 2 ) space = ^^k 2 AkS 2 W(k R). (3.6) 
IV. ^ FOR NON-GAUSSIAN FLUCTUATIONS 

In this section, we would like to recalculate the temperature fluctuation for the same 
power spectrum (2.3) but now for the NG configuration (2.4). Recall that in the G case the 
best one can do was to evaluate rms of various quanties of interest. For a NG case this is 
no longer necessary. For NG perturbations it is possible to get an explicit expression for 



a temperature fluctuation as a function of position x in a direction q. Using the fact that 
e i<t>{k) _ i ft j g trivial to do the angular sum dVt k in Eq. (2.6). 



AT 



q, x, r ) = —2n H^—^j J dkp(k) 1 ^ 2 (j (k\x + T q\) — monopole — dipole) . (4.1) 



T 

Using a thin Ak approximation and Eq. (3.2) with A = kr Q , A = q, B = kr, and B — x, to 
identify the monopole and the dipole terms, we get 

AT M,^) = -^f^W ( ) • (4-2) 



T ' ° 7 °V2W ° \+Sj 1 (k r )j 1 (k r)p 1 (q-x) 

where now q = \r q + rx\. It should be noted that this is the result for the specific NG 
configuration (2.4), and other NG configurations will give different results. 

V. ^-(R) FOR NON-GAUSSIAN FLUCTUATIONS 

We would now like to derive a simple expression for a mass fluctuation at some scale 
R that is appropriate for NG fluctuations. Again, for a NG perturbation it is no longer 
necessary to evaluate rms of quantities of interest. By definition (2.1), for (2.4) 

8p,^ / L N 3 



-(£) = (— ) / k 2 dkdn k \5 k \ 



p 

The angular integral dQ k is simple to evaluate. 



e i<f>(k) 



e lk - x . (5.1) 



Defining j{f(R) as an excess mass within a sphere of radius R centered on a point where the 
density contrast is maximum, we get 

S -^(R) =4tt(^)7 k 2 dk\5 k \y/W(kR), (5.3) 

where W{kR) is the familiar top down window function (3.5). Evaluating this for power 
spectrum (2.3) we get 



W {R) = ^(^yk 2 Ak5 ^W(k R). (5.4) 
VI. COMPARISON OF TWO TYPES OF PERTURBATIONS 

We will now use the simplest method to compare the magnitudes of two ^f- as predicted 
by each model. We will first write each AZ NG ( G ) j n terms of the corresponding ^-(R) NG ^ . 
Then we will compare the magnitudes of ^f NG and by setting ^-{R) NG = 6 -^{R) G '■ 
Using (3.3) with (3.6), we get for the G case 



(6.1) 



{q,A2) = \w{k R)~^Hlkf (^(R)) 

\ / i 



with q — r \qi — q 2 \- For the NG case using (4.2) with (5.4), we get 



s -±t ~ 3o(k T )j (k r) \ 
+3ji{k T )j l (k r)p l (q ■ x) J 

with q = \r q + rx\. The graph of for G case is shown in Figure 1. The graphs of 
for the NG case are shown in Figures 2(a)-2(f) for various values of k a r. It is interesting to 
note that for the NG case the temperature anisotropy depends on position in space. This 
is because the information about the positions of mass perturbations is not lost in the 
calculation for a NG case, whereas in the corresponding calculation for the G case they are 
lost by incorporating the rms methods. 

It should be noted that we are comparing the magnitudes of rms of autocorrelation 
function for ^ with the local temperature fluctuation; therefore, the angular parameter 9 
that appears in the figures for both G and NG have different meanings. 

In graphing ^jr-'s, a common normalization was used for both G and NG cases. From 
the figures we can conclude that, for the models described by (2.3) and (2.4), the predicted 
Sachs- Wolfe effect for the NG case is about 3 • 10" 3 smaller than for the G case. Actually 



AT 



{^) = -\w{k R)-^Hlk?^ (R) 



this conclusion is true as long as \r q — rx\ ^> k~ x . In other words, the position of spherical 
overdensity (void) does not intersect the locus of last scattering surface in relation to the 
observer at x. 

VII. DISCUSSION 

Before getting into the discussion of the generality of the result, some physical explana- 
tions are in order. First, the difference between G and the NG configuration (2.4) on ^-(x) 
is obvious. For the NG case, there are regions in space where the amplitudes of different 
wavelength fluctuations are adding coherently. In these regions the magnitudes of are 
larger than the rms of which is a relevant quantity for the G case. With this under- 

standing, the conclusion that the Sachs- Wolfe effect for the NG configuration (2.4) is smaller 
than the Sachs- Wolfe effect for the G configuration may seems contrary to Eq. (2.5). After 
all, if matter is more 'clumpy' in the NG configuration (2.4), than the maximum difference 
in the Newtonian potential should be larger for the NG case. And therefore, for the NG 
configuration should also be larger. The loophole in this argument lies in two facts. First, 
(2.5) needs to be normalized because of the quantity 8 (see (2.3), (3.3), (4.2)). Second, in 
calculating the Sachs- Wolfe effect I = (monopole) term of ¥f- is subtracted out because 
it is unobservable. For the NG configuration (2.4) under consideration, it is easy to see 
that regions in space where ^-(q, x) is much larger than the rms of ^-(q, x) are spherically 
symmetric. Therefore, when the / = term is subtracted out the enhancement of going from 
a linear sum to the 'random walk' sum largely disappears. 

We would now like to give a proof that the NG configuration (2.4) among all possible NG 
configurations gives the smallest Sachs- Wolfe effect for the power spectrum (2.3). As stated 
in the introduction, even for a fixed power spectrum, the set of all NG field configurations is 
an infinite parameter space. For the argument we choose to parametrize a NG configuration 
in a following way. A general phase of a fluctuation 5k can be expanded in terms of spherical 
harmonics 



e im = Y,YU'k)GU\k\)- 



(7.1) 



l,m 



For the power spectrum (2.3) we are considering, \k\ is constant, so G^ m {k) are just constants. 

. , i 1 2 i ■ I 

Since e l(p is a unit norm, it implies that J2 \Gi,m\ = 47r or \Gi jTn \ < v47r. In this general NG 
field configuration, it is easy to show that for fixed power spectrum the new mass fluctuation 
is (Eqs. (5.1)-(5.4)) 



SM 



(R) 



NG 



G , 5M_ 
4¥ M 



(R) 



old 



NG 



(7.2) 



By old we mean the previously discussed case e l<t> = 1. We are being little cavalier here. 
Even though the resulting density configuration ^-(x) may not be spherically symmetric, we 
are defining ^-(R) as an excess mass within a sphere of radius R centered on a point where 
the density contrast is maximum. 

The change in the Newtonian potential is more complicated. Inserting (7.1) into (2.6) 
and after dVLk we get 



5<t>{x + r q) = -6nH 2 (^) Ak8 ^i% m (t)3i(ko<;)G l , r , 

V27r/ l,m 

with <f = x + r q, and x) = h[6<f>(s) — monopole] with (see (2.5)) 



(7.3) 



monopole term 



5(j)(x)j (k r). 



(7.4) 



Denoting quantities with the S Q and G^ m dependence extracted by a hat (e.g., A) we get 
from (7.3), (7.4), and (2.5) 



-jr(Q,x) = jY,Gi, m m, m (3 - S<f> ltm (x)j (k r) 



(7.5) 



And from (5.4) and (7.2) we get ^ = 5 G 00 ^-. Eliminating 5 Q , we need to find the 



extremum of 



l,m 

2 

subject to the constraint X) |G; )m | = 47r. With the method of lagrange multiplier, we get 
X) |G7,m| 2 = G 2 00 i or e*^ = G 00 Y 00 = ±1. It is easy to check that the extremum is a minimum 
because another configuration, e.g., G OQ 0, gives a larger value. It should be noted that 
the variational method alone does not tell us that the field configuration e l<t> = ±1 gives a 
Sachs- Wolfe effect that is smaller than for the G configuration. A remark on the measure 
of generality of the result is possible, given by Eq. (7.6) is a smooth linear function of 
Gi )Tn except for G OQ . And therefore, combined with the restriction of \Gi tm \ 2 = An, the 
conclusion will be unchanged as long as G 0Q ~ ±v^4vr. 

The result for a more realistic power spectrum p{k) ~ k n will be discussed elsewhere [42]. 
The decrease in the Sachs- Wolfe effect persists but to a different degree. 

In light of the result that there are NG field configurations that can reduce the Sachs- 
Wolfe effect, we could turn the analysis on its head. By properly incorporting void evolution, 
we could have addressed the question of the largest possible void compatible with COBE 
given a NG density fluctuations (here e 1 ^ ~ —1). This would be the NG counterpart of the 
work by Blumenthal et al. [43]. 

And finally, even though in our analysis we have assumed that perturbations and the scale 
factor of the universe evolved as that of flat FRW universe, the decrease in the Sachs- Wolfe 
effect will obviously persist in other cosmologies but to different degrees. 

IIX. CONCLUSION 

In this paper, we have taken a power spectrum, which was essentially a delta function, and 
have shown that there is a set of non-Gaussian field configurations for which the Sachs- Wolfe 
effect is smaller than the corresponding effect for the Gaussian configuration case. This set of 
non-Gaussian perturbations is characterized by having a large monopole term , G ao ~ ±\^Atv, 
in the Y^ m decomposition of the phases, i.e., e 1 ^ = J2i,mYi,m(k)Gi t7n (\k\) xs ±1. 



The conflict between the recent COBE measurements of CMBR and some HDM Gaussian 
models can potentially be salvaged by the combination of a non-Gaussian fluctuation model 
and HDM. 



APPENDIX 

We would like to prove an identity. 



U\AA -BB\)= £(2Z + l^A^B^A ■ B) 
i 



consider 



da 



J e KAA-BByc dQc = J^t l (2l + l)j l (\AA-BB\) Pl ((AA-BB)-C)dQ c 
= 4ir Jo (\AA-BB\) 

= /E ^( 2/ + l )3i( A M A ■ C)] [(-if (21' + l)ji>(B) Pl ,(B ■ C) 

using J Pl (A ■ 6)p v (B ■ C)dn c = ^jPi(A ■ B)8 W 

= 47rY,(2l + l)3i(A)ji(B) Pl (A-B). Q.E.D 
i 
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FIGURE CAPTIONS 

1. Graph of the square root of the temperature autocorrelation function | ^ AT X e °) AT(e<,+e) ^ 
Eq. (6.1), as function of 9. H' 1 = 3000MPC and lTik~ x = 40MPC was used to imi- 
tate the peak in the power spectrum for a HDM. The normalization is arbitrary. 

2. Graph of the temperature fluctuation ^-(q, x), Eq. (6.2), as function of 9 = cos~ l (q- x) 
for various values of r = \x\. Figures 2(a)-2(f) correspond to k a r = [0, 2.0, 4.0, 6.0, 8.0, 1( 
respectively. The values of H Q , k Q , and the overall normalization are the same as those 
of Fig. 1. 
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